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Abstract. We have recently presented a collision- model-based framework to 
approach non-Markovian quantum dynamics [Ciccarello F Palma G M and Giovannetti 
V 2013 Phys. Rev. A 87 040103(R)]. As a distinctive feature, memory is introduced in 
a dynamical way by adding extra inter-ancillary collisions to a standard (memoryless) 
collision model. Here, we focus on the case where such intra-bath collisions arc 
described by incoherent partial swap operations. After briefly reviewing the model, 
we show how to include temperature as an additional parameter by relaxing the 
assumption that each bath ancilla is initially in a pure state. We also calculate explicitly 
the dynamical map entailed by the master equation in the paradigmatic instance of a 
Jayncs-Cummings system-ancilla coupling. 
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1. Introduction 

The theoretical settling of open system dynamics ranks among the hottest problem in 
modern quantum mechanics [USE]. By definition, an open system is one in contact with 
an external environment and thus the dynamical map describing its time evolution is 
non-unitary. On a rather general basis, a physically grounded dynamical map is required 
to be completely positive and trace preserving (CPT). In the case that the environment 
is memoryless, a Markovian dynamics occurs which, demonstrably, is always governed 
by a so called Lindblad-type master equation (ME) [H [21 [3]. Any Lindblad ME in 
turn corresponds to an unconditionally CPT dynamical. In the general case, though, 
environments are intrinsically non-Markovian (NM) and in several known scenarios a 
Markovian-based approach is fully inadequate [I]. Unlike the Markovian case, however, 
a general systematic framework to tackle NM processes has not been developed to date 
[5]. Rather, a number of variegated approaches have been proposed. Usually, the major 
drawback of these is that the corresponding MEs can violate the CPT condition in 
some regimes [61 [7J Ej- Typically, this stems from the phenomenological assumptions 
and approximations that one somehow needs to make in order to derive a ME j9]. A 
further recurrent drawback is that some MEs, such as those in Refs. [TOl [TTj . focus on 
a regime which is too weak a perturbation of pure Markovianity. As a result, their 
capability to capture genuinely NM features can be severely limited [T2] . 

Recently, we have proposed an innovative collision-model-based approach to 
tackle non-Markovian dynamics [13]. Collision models (CMs) were first envisaged by 
Boltzmann as an interesting way to describe the irreversible dynamics of a system in 
contact with a large environment, or bath, in terms of successive "collisions" with its 
small subparts (ancillas). In more recent years, mostly, they were successfully applied 
to study quantum open systems [HJ HSl [161 HI EH]- I n a standard CM, the bath is 
modeled as a large collection of non-interacting identical ancillas. S (the system under 
study) "collides" with each of these one at a time. It can be shown that such a process 
gives rise to an irreversible dynamics for S exactly described by a Lindblad-type ME 
[TT1 ITS] . This is because S is not allowed to interact more than once with a given 
ancilla. The bath thus cannot keep track of the system's past history. A major reason 
why CMs are attractive is that they are suited to work out exact MEs basically without 
any approximations (hence ruling out the possibility to violate the CPT condition). 
Only the passage to the continuous limit is needed [131 El UHl [19] . This is in contrast to 
standard microscopic system- reservoir models [TJ, where even to derive Markovian MEs 
drastic assumptions (such as the Born-Markov approximation) are in fact unavoidable. 
Very recently, Rybar et al. introduced a NM CM able to simulate any indivisible single- 
qubit channel [19]. In their framework, memory is introduced by simply allowing for 
correlated initial bath states (aside from this, the model is identical to memoryless 
CMs). 

In our model [13], instead, the initial reservoir state is fully uncorrelated (just like 
in a standard Markovian CM) but ancillas can undergo pairwise collisions. This enables 
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transmission of quantum information across the bath in a way that the initial state of a 
given ancilla before colliding with the system is affected by the S past history. Moreover, 
this mechanism endows the reservoir with memory in a dynamical way: the ability to 
remember is associated with physical parameters entering the system-bath interaction, 
which is what one would intuitively expect. In Ref. [13], we have proposed to describe 
each inter-ancillary collision as a partial swapping operation. This is indeed a natural 
choice to account for intra-bath information transfer. Moreover, it is defined in terms 
of the swap unitary operator, which allows to considerably simplify calculations. A 
possible choice is what we call an incoherent swap, i.e, a map whose Kraus operators 
are proportional to the identity and swap operators, respectively [13] . Using this, which 
is non-unitary but CPT, we have demonstrated that one can work out an exact ME 
interpolating between a fully Markovian regime and a strongly NM one. Importantly, 
such ME entails a dynamics which is unconditionally CPT [13]. Furthermore, it does not 
depend either on the system-ancilla coupling form or the dimensionality of the particles. 

In this short paper, we consider the aforementioned incoherent swap CM and 
present some related developments with a twofold goal. First, our theory in Ref. [T3] 
was restricted to the case that the initial bath state is pure. Yet, any open system 
model is expected to feature temperature among its parameters. This is usually taken 
into account by considering the reservoir initially in a thermal, thus in general mixed, 
state. In this respect, we show that assuming an initial pure state is not restrictive 
provided that one redefines the ancillas as initially entangled bipartite systems where 
only one of the two can interact with S (an approach that in the theory of quantum 
channels [20] corresponds to the Stinespring dialation mechanism). Our second aim is to 
illustrate the explicit form of the solution of our ME in a paradigmatic case. We choose 
the situation where the system-ancilla coupling is of the Jaynes-Cummings (JC) form. 
This, indeed, routinely occurs in quantum optics and is often used as an illustrative 
instance of various proposed approaches to NM dynamics. 

This paper is structured as follows. In Section [21 we briefly review the theory 
leading to the discussed ME for incoherent swap intra-bath collisions. In Section [31 
we show how to include temperature in our framework. In Section HI we consider a 
system-ancilla interaction having the form of a JC coupling and explicitly work out the 
dynamical map entailed by our ME. Finally, in Section EJ we draw our conclusions. 

2. Review of the collision model for incoherent partial swap 

Our model comprises a system S, initially in state po; an d a bath consisting of a collection 
of ancillas labeled with i — 1,2, .... The initial overall state is <To = PoPbo, where pbo is 
the bath initial state given by 



where |0) B = |0) 1 |0) 2 ---. As mentioned, in our model inter-ancillary (AA) collisions are 
interspersed with system-ancilla (SA) ones. In the beginning, S collides with ancilla 1. 




(1) 
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In standard (Markovian) CMs, 5-2 collision would then follow, then 5-3 and so on. 
This way, each ancilla would still be in the initial state |0)(0| before colliding with 5, 
thus fully "unaware" of previous collisions. In contrast, in our model we assume that an 
extra AA collision between 1 and 2 occurs after 5-1 but before 5-2. Thereby, prior to 
its interaction with S, ancilla 2 will now be in a perturbed state in which information 
over past history of S is imprinted. The process proceeds by mere iteration: once S-2 
collision is over, a 2-3 interaction follows, then 5-3, 3-4 etc. 

Each collision, either SA or AA, is described by a CPT quantum map affecting the 
degrees of freedom of the two involved particles. Specifically, the map for an SA collision 
involving the ith ancilla is defined as the unitary coupling UsiO~ = e~ lHsitc ae tHsitc , where 
t c and Hsi are respectively the collision time and interaction Hamiltonian (we set h — 1 
throughout). Instead, the AA collision involving the ancilla % to the jth one is defined 
in terms of the CPT map 

Sija = (l - p s )a+p s S ij aS ij (2) 

where p s is the swap probability while S^ is the well-known swap operator [2Tj 
exchanging the states of i and j as S^ ly^JV^j = I^Jv 9 )^ The steps of our discrete 
process are defined in a such way that the nth step terminates when both (n-l)-n and 
S-n collisions are over (but step n = 1 ending after 5-1 collision). Hence, at the nth 
step the overall state is given by a n = (Us n °<Sn,n-i° ■ ■ ■ °Us2 S2,i°Usi) o" (henceforth, 
the symbol "o" which describes super-operator composition will be omitted). Let 
p n = Tr B (a n ) be the nth -step state of 5 (the partial trace is over all the ancillas). Using 
the properties of the swap operator S iji+ ie~ lHsitc S i:i+ i = e - tH s,n-itc an( j g.. |0) B = |0) B , 
one can easily show [13] that p n can be expanded in terms of all previous-step states 
of 5 only. In the continuous limit, where due to t c ~ and n ^> 1 the step number is 
replaced by the continuos time t = nt c , such expansion for p n can be shown to give rise 
to a corresponding integro-differential ME. This reads 

Tt -TJ Q dte £ {t)-^-^ + e —po (3) 

with the CPT map S(t) defined as 

£(*)p=£nHe-*iW|0) n p n v \ e -4&Bnt\ 0) \i } (4) 



where each Kraus operator n (y\ e~ tHsnt \0) n is independent of n since so is the form 
of Hsn and the ancillas are all identical ({|f)} is a basis for the single-ancilla Hilbert 
space). Map S(t) is in general strongly NM since, evidently, it effectively describes 
the continuous coherent interaction of 5 with a single ancilla. Indeed, our model 
interpolates between two extreme regimes depending on the value of T. When T = 0, 
which can be shown to correspond to p s — 1 in the discrete model, swapping is perfect 
and as is immediate to see from Eq. ([3]) the solution reads p(t) = S(t)p . In such 
a case, then, Eq. (j4j) coincides with the process dynamical map A(t) [this is defined 
through p{t) = A(t)po\. In the opposite limit r — > oo, instead, it can be proven [13] 
that A(t) = e £ ^ 1 with Eq. ([3]) reducing to the Lindblad form p = £(0)p, where £(0) 
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is a Lindlabian superoperator. Such limiting case indeed corresponds in the discrete 
picture to p s = 0, i.e., AA collisions do not occur at all [cf. Eq. (J2J)] and a standard fully 
Markovian CM is retrieved. 

The general expression for A(£) can be worked out as follows. Clearly, A(t) obeys 
Eq. ([3]) under the formal replacement p— >A. By taking the Laplace transform (LT) of 
such equation, this is easily solved as [13] 



where A(s) and 8 (s) are the LTs of A(t) and 8{t), respectively. Expanding Eq. (J5]) in 

~ r ~ i k 

powers of V gives A(s) = Y,T=i £( S + T) I 1 * 1-1 , whose inverse LT is 



A(t)=£" 1 [A( a )](t) =Er H C-\i h (8+T)](t) . (6) 

fc=l 

Basic properties of LT allow to immediately calculate the inverse LT on the right-hand 
side as 

£- 1 [^+r)] = e- r * f dtl f l dh ... f^dt^ 8{t k ^8{t k _ 2 -t k ^) ■■■S(t-t 1 ). (7) 
Jo Jo Jo 

We have thus expressed A(t) as a weighted series of multiple auto-convolutions of the 
CPT map 8(t). Being a composition of CPT maps, each convolution [the integrand in 
Eq. ([7])] is CPT itself. Also, it is multiplied by a positive coefficient [cf. Eqs. <Q and 
(jJJ)], which yields complete positivity of map A(i). Moreover, it is easily checked [13] 
that regardless of T and S(t) Tr [A(£)po] = 1- Map A(t) is thus always CPT. 

3. Including temperature 

As mentioned in the Introduction, the natural way to account for temperature T is to 
replace the bath initial state Eq. (CQ) with a product of single-ancilla thermal states as 

where (3 = 1/(KT) {K is the Boltzmann constant) and Hi is the single-ancilla free 
Hamiltonian (whose operator form is independent of i). For d- dimensional ancillas, 
let (k = 0, ...,d— 1) be the orthonormal set of eigenstates of Hi such that 

H i \n) i = e K Then, the ith-ancilla thermal state p^o is i n general a mixture of 

according to 

dr-l e -f3s K 

Pio = Yl \ K )L K \ > ( 9 ) 

k=0 z 

where Z = J2t~lo e ~^ £K ^ s the partition function. Clearly, for (3 — > oo, namely zero 
temperature, state Eq. (jSJ) reduces to a product of pure states such as Eq. ([1]). For 
non-zero temperatures, the ancillary initial state can be "purified" [21] as follows. In 
general, given a single-particle mixture describing the state of a particle A, one can 
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define a further fictitious particle B in a way that A-B are in a pure global state AB 
and = Tr s (| , 0) j4B ('i/'|). Following this scheme, we thus replace in our model each 
d-dimensional ancilla with a pair of identical d- dimensional ancillas such that only one 
of the two can interact with S (d is an arbitrary integer). We now use indexes i\ 
and «2 to label the ancillas really interacting with 5* and the corresponding fictitious 
ones, respectively. Each ixih. ancilla is initially in state p il0 as given by Eq. (Q. 
Correspondingly, the joint system comprising ii and the associated auxiliary ancilla 
22 is taken to be initially in the pure state 

l^o) il!2 = Ej— l fi )J K ) i2 > (io) 

which evidently fulfills pj 1 o = Trj 2 (|^ )(' l /'o|)- 

In the derivation of Eq. fl3]) no assumption on either the dimensionality of the 
ancillas or the features of H$i is necessary. Therefore, map S(t) now reads 

£(t)p=J2 M e ~ iAsilt l^o) P (W 2 | e- i6s ^ \^ )y , (11) 

where {l^i)}^ ({|^2)}i 2 ) I s a basis for i\ (i 2 ). Using that Hsh does not affect we end 
up with 

£{t)P = T, {E H e- <JW l«> P (Wile-^ |«)) f } , (12) 

where our notation emphasizes that the resulting map is a convex combination of maps 
defined in terms of initial pure single- ancilla states {]/«)}■ We conclude that ME Eq. ([3]) 
holds at finite temperature as well provided that map S(t) in Eq. ( fT2l) is considered. 
Note that Eq. (I12p also shows that S(t) can now be regarded as the map describing the 
reduced unitary dynamics associated with the time-evolution operator e~ tHsi i l and the 
initial S-i\ state given by the tensor product of p with Eq. (jUJ). Hence, it is still true 
that for r = S effectively behaves as it is continuously interacting with a single ancilla. 

4. Dynamical map for Jaynes-Cummings coupling 

We assume that S is a qubit [21], i.e., a two-level system, and the SA coupling has 
the form of an XY isotropic interaction corresponding to a JC model when S and 
A are on resonance [22]- Each ancilla is modeled as a bosonic mode initially in the 
vacuum state, but because of conservation of the total number of excitations it behaves 
as an effective qubit as well. Let then {|0) S /^,|1) S /^} be a basis for S (the zth ancilla). 
The SA Hamiltonian reads Hsi — ^ ftT+5i_+H.c), where Q is a coupling rate while 
a + = = |1) 5 (0| and S i+ = = |l),j(0| are the usual spin-1/2 ladder operators. The 
most general initial state of S reads po = (l— p) |0) 5 (0|+p |l) s (l|+(r |0) 5 (l|+H.c), where 
<p< 1 and |r| 2 <p(l— p). An amplitude damping channel (ADC) [21] transforms po 
into 

A(v)po=[l ~ V 2 P] |0) 5 (0| + rfp |l) s (l| + [nr \0) s (l\ + H.c] , (13) 
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where <rj < 1. Through a standard calculation, when the ancillas are initially all in |0) 
the map in Eq. (T4J) is given by S(t) = A[cos(Qt)]. ADCs fulfill the composition property 
A(r}i)A(r}2)=A(rjir]2)- Using this in combination with Eqs. (EJ) and ([7]) and introducing 
the rescaled time r = fit, the dynamical map can be expressed in matrix form as 

' /?i(r)r \ 

/9i(r)r* /? 2 (r)p J 



A(r)p -- 



(14) 



with 

AW 



oo 

-rr ^ pjfc-l 
fe=l 



dr± [ dr2 



Tk-2 



dr k „i cos(tv-i)^ cos(r fc __ 2 -r fc _i) £ •••cos(r-ri)^, (15) 



where f = T/f2. Hence, functions /?i(r) and /^(r) fully specify the dynamical map. We 
define the LT of a cosine power as q(s) = £[cos(t)*]. Given that Eq. (JIB]) is formally 
analogous to Eq. (jHJ) with due replacements, in line with Eq. (jHJ) the LT of (3i(t) is 
evidently given by 

Ci(s+T) 



Pi(s)=£\p t (T)]: 



(16) 



1-Td e (s+T) 

Now, by using Ci(s) = s/(s 2 +l) along with C2(s) = (s 2 +2)/[s(s 2 +4)] and anti-transforming 
we find 



Tt 

■e 2 



/3 2 (r)=5> ie 



rsinh(iVr 2 -4r 



VP 



cosh 



P-4r 



;i7) 



i=i 



where 



«i 



r - Vr 3 + 35 + 9r 



12 



3VT 3 + 35 + 9r 



02 = 0^ = - 

3 g 



v/3 + 



r 3 + 35 + 9r 



i + 1V3) fr 2 - 12 



Vr 3 + 35 + 9r 



4r 



A x -. 



2r«i + a\ + r 2 + 2 
cc 2 + |«2| 2 — 2aiRe(a2) 



Ao = Al 



i (2r«2 + ot 



2(ai-a2)Im(a 2 ) 



with 5 = V 6r 4 — 39r 2 + 192. We have checked that /3i(r) and /3 2 (r), as given by 
Eqs. (ED and ([181), fulfill the conditions < /3 2 (r) < 1 and /?i(r) 2 < /3 2 (r) regardless 
of T and r. This confirms that map Eq. (fl4|) is unconditionally CPT in accordance with 
our general theory. 



5. Conclusions 



In conclusion, we have briefly reviewed a CM-based framework to tackle non-Markovian 
dynamics. The model differs from standard Markovian CMs in that extra inter-ancillary 
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collisions are added between next system-ancilla interactions. Here, we have focussed on 
the case where each ancilla-ancilla collision is described by an incoherent partial swap. 
The corresponding dynamics is exactly described by a unconditionally CPT ME. We 
have complemented the related theory with two further developments. First, we have 
shown that the case where the bath is initially in a thermal state, instead of a pure one, 
is in fact already encompassed in the theory provided that one suitably redefines each 
ancilla as a bipartite system (initially entangled). Next, we have illustrated for the first 
time a paradigmatic instance where the solution of the ME can be calculated explicitly 
in a compact form. We have chosen the case where the system-ancilla interaction is 
described by a JC coupling. 
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